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ABSTRACT 



We investigate the self-consistent electrodynamic structure of a particle accelerator in the 
Crab pulsar magnetosphere on the two-dimensional poloidal plane, solving the Poisson equation 
for the electrostatic potential together with the Boltzmann equations for electrons, positrons and 
gamma-rays. If the trans-field thickness of the gap is thin, the created current density becomes 
1/-^ \ sub-Goldreich- Julian, giving the traditional outer-gap solution but with negligible gamma-ray 

' luminosity. As the thickness increases, the created current increases to become super-Goldreich- 

. Julian, giving a new gap solution with substantially screened acceleration electric field in the 

\^ I inner part. In this case, the gap extends towards the neutron star with a small-amplitude 

■ positive acceleration field, extracting ions from the stellar surface as a space-charge-limited flow. 

The acceleration field is highly unscreened in the outer magnetosphere, resulting in a gamma-ray 
Oh' spectral shape which is consistent with the observations. 

Y^, I Subject headings: gamma-rays: observations - gamma-rays: theory - magnetic fields - methods: 

' numerical - pulsars: individual(Crab) 

> 
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' 1. Introduction 

The Energetic Gamma Ray Experiment Telescope (EGRET) aboard the Compton Gamma Ray Ob- 
servatory has detected pulsed signals from at least six rotation-powered pulsars (e.g., for the Grab pulsar, 
Nolan et al. 1993, Fierro et al. 1998). Since interpreting 7-rays should be less ambiguous compared with 
reprocessed, non-thermal X-rays, the 7-ray pulsations observed from these objects are particularly important 
as a direct signature of basic non-thermal processes in pulsar magnetospheres, and potentially should help 
to discriminate among different emission models. 

The pulsar magnetosphere can be divided into two zones (fig. 1): The closed zone filled with a dense 
plasma corotating with the star, and the open zone in which plasma flows along the open field lines to escape 
through the light cylinder. The last-open field lines form the border of the open magnetic field line bundle. 
In all the pulsar emission models, particle acceleration and the resultant photon emissions take place within 
this open zone. 
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On the spinning neutron star surface, an electro-motive force, EMF ?a fi^B+rJ/c^ « lO"'^^'^ V, is exerted 
from the magnetic pole to the rim of the polar cap. In this paper, we assume that both the spin and magnetic 
axes reside in the same hemisphere; that is, • /ii > 0, where represents the rotation vector, and /x the 
stellar magnetic moment vector. This strong EMF causes the magnetospheric currents that flow outwards in 

the lower latitudes and inwards near the magnetic axis (left panel in fig. 2). The return current is formed at 
large-distances where Poynting flux is converted into kinetic energy of particles or dissipated (Shibata 1997). 

Attempts to model the particle accelerator have traditionally concentrated on two scenarios: Polar-cap 
models with emission altitudes of ~ lO^cm to several neutron star radii over a pulsar polar cap surface 
(Harding, Tademaru, & Esposito 1978; Daugherty & Harding 1982, 1996; Dermer & Sturner 1994; Sturner, 
Dcrmer, & Michel 1995), and outer-gap models with acceleration occurring in the open zone located near 
the light cylinder (Cheng, Ho, & Ruderman 1986a, b, hereafter CHR86a,b; Chiang & Romani 1992, 1994; 
Romani and Yadigaroglu 1995). Both models predict that electrons and positrons are accelerated in a charge 
depletion region, a potential gap, by the electric field along the magnetic field lines to radiate high-energy 
7-rays via the curvature and inverse-Compton (IC) processes. 

In the outer magnetosphere picture of Romani (1996), he estimated the evolution of high-energy flux 
efficiencies and beaming fractions to discuss the detection statistics, by considering how pair creation on 
thermal surface flux can limit the acceleration zones. Subsequently, Cheng, Ruderman and Zhang (2000, 
hereafter CRZOO) developed a three-dimensional outer magnetospheric gap model, self-consistently limiting 
the gap size by pair creation from collisions of thermal photons from the polar cap that is heated by the 
bombardment of gap-accelerated charged particles. The outer gap models of these two groups have been 
successful in explaining the observed light curves, particularly in reproducing the wide separation of the two 
peaks commonly observed from 7-ray pulsars (Kanbach 1999; Thompson 2001), without invoking a very 
small inclination angle. In these outer gap models, they consider that the gap extends from the null surface, 
where the Goldreich- Julian charge density vanishes, to the light cylinder, beyond which the velocity of a 
co-rotating plasma would exceed the velocity of light, adopting the vacuum solution of the Poisson equation 
for the electrostatic potential (CHR86a). 

However, it was analytically demonstrated by Hirotani, Harding, and Shibata (2003, HHS03) that an 
active gap, which must be non- vacuum, possesses a qualitatively different properties from the vacuum solution 
discussed in traditional outer-gap models. For example, the gap inner boundary shifts towards the star as 
the created current increases and at last touch the star if the created current exceeds the Goldreich-Julian 
(GJ) value at the surface. Therefore, to understand the particle accelerator, which extends from the stellar 
surface to the outer magnetosphere, we have to merge the outer-gap and polar-cap models, which have been 
separately considered so far. 

In traditional polar-cap models, the energetics and pair cascade spectrum have had success in reproduc- 
ing the observations. However, the predicted beam size of radiation emitted near the stellar surface is too 
small to produce the wide pulse profiles that are observed from high-energy pulsars. Seeking the possibihty 
of a wide hollow cone of high-energy radiation due to the flaring of field lines, Arons (1983) first examined 
the particle acceleration at the high altitudes along the last open field line. This type of accelerator, or 
the slot gap, forms because the pair formation front (PFF), which screens the accelerating electric field, 
i^ll , in a width comparable to the neutron star radius, occurs at increasingly higher altitude as the mag- 
netic colatitude approaches the edge of the open field region (Arons & Scharlemann 1979). Muslimov and 
Harding (2003, hereafter MH03) extended this argument by including two new features: acceleration due 
to space-time dragging, and the additional decrease of at the edge of the gap due to the narrowness of 
the slot gap. Moreover, Muslimov and Harding (2004a, b, hereafter MH04a,b) matched the high-altitude slot 
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Fig. 1. — Schematic figure (sideview) of tlie two representative accelerator models. The small filled circle 
represents the neutron star 




Fig. 2. — Schematic picture of electric current in the pulsar magnetosphere. Left: Global electric current 
due to the EMF exerted on the spinning neutron star surface when fl ■ fi > 0. Right: Current (downward 
arrows) derived in the inncr-slot-gap (shaded region). 
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gap solution for to the solution obtained at lower altitudes (MHOS), and found that the residual is 
small and constant, but still large enough at all altitudes to maintain the energies of electrons, which are 
extracted from the star, above 5 TeV. 

It is noteworthy that the polar-slot gap model proposed by MH04a,b is an extension of the polar-cap 
model into the outer magnetosphere, assuming that the plasma flowing in the gap consists of only one sign of 
charges. This assumption is sclf-consistently satisfied in their model, because pair creation in the extended 
slot gap occurs at a reduced rate and the pair cascade due to inward-migrating particles docs not take place. 
In the polar-slot gap model, the completely charge-separated, space-charge-limited flow (SCLF) leads to a 
negative E^^ for O • /x > 0. However, we should notice here that the electric current induced by the negative 
(right panel of figure 2) contradicts with the global current patterns (left panel of figure 2), which is 
derived by the EMF exerted on the spinning neutron-star surface, if the gap is located near the last-open 
field line. (Note that the return current sheet is not assumed on the last-open field line in the slot gap 
model.) 

On these grounds, we are motivated by the need to contrive an accelerator model that predicts a con- 
sistent current direction with the global requirement. To this aim, it is straightforward to extend recent 
outer-gap models, which predict opposite i^y to polar-cap models, into the inner magnetosphere. Extend- 
ing the one-dimensional analysis along the field lines in several outer-gap models (Hirotani and Shibata 
1999a, b, c; HHS03), Takata, Shibata, and Hirotani (2004, hereafter TSH04) and Takata et al. (2006, here- 
after TSHC06) solved the Poisson equation for the electrostatic potential on the two-dimensional poloidal 
plane, and revealed that the gap inner boundary is located inside of the null surface owing to the pair creation 
within the gap, assuming that the particle motion immediately saturates in the balance between electric and 
radiation-reaction forces. 

In the present paper, we extend TSH04 and TSHC06 by solving the particle energy distribution explicitly, 

and by considering a super-GJ current solution with ion emission from the neutron star surface. In § 2, we 
formulate the basic equations and boundary conditions. We then apply it to the Crab pulsar in § 3, and 
compare the solution with MH04 in § 4. 

2. Gap Electrodynamics 

In this section, we formulate the basic equations to describe the particle accelerator, extending the 
method first proposed by Beskin et al. (1992) for black-hole magnetospheres. 

2.1. Background Geometry 

Around a rotating neutron star with angular frequency SI, mass M and moment of inertia I, the 
background space-time geometry is given by (Lense & Thirring 1918) 

ds'^ = gttdt^ + 2gt^dtdip + grrdr^ + geedO'^ + g^^pdip^, (1) 

where 

gtt = (l - ^) c2, gt^ = ac^ sin^ (9, (2) 

grr = - (l - — ) , 9ee = -r^, g<pyp = -r^ sin^ 6*; (3) 
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Tg = 2GM/(? indicates the Schwarzschild radius, and a = IQ/{Mc) parameterizes the stehar angular 
momentum; second and higher order terms in the expansion of a/r^ are neglected. At radial coordinate r, 
the inertial frame is dragged at angular frequency 

where r* represents the stellar radius, I45 = //lO^^ erg cm^, and tq = r/lOkm. 



2.2. Poisson Equation for Electrostatic Potential 

The first kind equation we have to consider is the Poisson equation for the electrostatic potential, which 
is given by the Gauss's law as 



(5) 



where V denotes the covariant derivative, the Greek indices run over t, r, 6, ip; = \J grrQeeP^ = cr^ sin 9 
and 



pI = gL - gttg<p<p = c2 ( 1 - ^ ) sin^ e. (6) 



r 



If there is an ion emission from the stellar surface into the magnetosphere, the total real charge density p is 
given by 

P = Pe + Pion, (7) 

where Pc denotes the sum of positronic and electronic charge densities, while pion does the ionic one. The 
six independent components of the field-strength tensor give the electromagnetic field observed by a distant 
static observer (not by the zero-angular-momentum observer) such that (Camenzind 1986a, b) 

Er = Frt, Ee = Ffft, E^ = F^t, (8) 

pr _ 9tt + gty^ p p0 _ 9tt + gty^ p „ Pw p /qn 

-D — ^ — i^Sip, -D — ^ — -^(pr, iJif = — ^ — ^ ^ 

where F^^^ = A^^^^ — and ^^,1/ denotes the vector potential partially differentiated with respect to 
a;". The strength of the poloidal component of the magnetic field is defined as 



Bp = c ■ . (W) 

gtt + gt<fiii' 

Assuming that the electromagnetic fields are unchanged in the corotating frame, we can introduce the 
non-corotational potential \1/ such that 

F^^t + OF^v = -d^^ir, e,<f- ^t), (11) 

where p, = t,r, 6, (p. If F^t + OFa,^ = holds for A = r,9, the angular frequency O of a magnetic field is 
conserved along the field line. On the neutron-star surface, we impose Fgt + flFg^ = (perfect conductor) 
to find that the surface is equi-potential, that is, dg'if = df'i' + ^d^'i/ = holds. However, in a particle 
acceleration region, F^t + ^F^^ deviates from and the magnetic field does not rigidly rotate (even though 
the deviation from the uniform rotation is small when the potential drop in the gap is much less than the 
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EMF exerted on the spinning neutron star surface). The deviation is expressed in terms of ^, which gives 
the strength of the acceleration electric field that is measured by a distant static observer as 



^\ = ^-E= ^{Fu + = I • (-V*), (12) 



where the Latin index i runs over spatial coordinates r, 6, (f, and an identity B^'Fr^ + B^Fq^ = is used. 
Substituting equation (11) into (5), we obtain the Poisson equation for the non-corotational potential 

^2 



( ^g'^'g^vd'^'if ] = 47r(p - pGj), (13) 



V-9 \ P- 

where the general relativistic Goldreich-Julian charge density is defined as 



2 -g''''gMn-u;)F^^ 

Pw 



(14) 



Using g-- = 1/grr, g^' = 1/gee, g'" = g'"- = 0, = -g^^/pl, 5^* = gt^/pl, and g^^ = -gu/pl, taking 
the limit r ^ Tg, and noting that dr{r^ —geeB^) — de{\/ —grrB'') gives the toroidal component of V x B, we 
find that equation (14) reduces to the ordinary, special-relativistic expression of the Goldreich-Julian charge 
density (Goldreich and Juhan 1969; Mestel 1971). 

Instead of {r,9,ip), we adopt in this paper the magnetic coordinates (s,0*,<^*) such that s denotes the 
distance along a magnetic field line, ^* and the magnetic colatitude and the magnetic azimuth of the 
point where the field line intersects the stellar surface. Defining that = corresponds to the magnetic 
axis and that (p* = to the plane on which both the rotation and the magnetic axes reside, we can compute 
spherical coordinate variables as follows: 

r(s0 ^)-r +r ?:i^2hl^ds' (15) 

e{s,e.,<p.) = e{Q,e^,^.)+ f ^^^Ai^^^^ds' , (i6) 

Jo B(s',e,ip-iit} 

Jo B(s',6',(^ - m) 

where 0(O,0,,(/3*) satisfies sin (?(0, (/?*) • cosi^H.sinai + cos 0(0, 0,, (/?*) • cosai = cos0*; a\ represents the 
angle between the rotation and magnetic axes. We can numerically compute the transformation matrix 
dx^ jdx^ and its inverse dx^ jdx^ from equations (15)-(17), where x^ =r, x^ = 6, x^ = ip, X'^ = s, x^ = 
and x^ = <^*. Substituting 

d _ ds d ^ 90* d ^ dip* d 
dx^ dx^ ds dx^ 90* 9a;' dipt, 

into equation (13), and utilizing d/dt = —fld/d(p, we obtain the following form of Poisson equation, which 
can be applied to arbitrary magnetic field configurations: 

(g-92 + g'-'^dl + g^'^'dl^ + 2/^*9.9,, + 25^*^*9^.9^. + 25^*^9^.9.) * 
- {A'ds + A«*9e. + A^'d^,) * = 47r(p - pcj), (19) 
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where (see Appendix for explicit expressions) 




fco / dx^ 



r,0 



dx^ 
dip 



(20) 



rfi 



and 
A}' = 



Or 



„2 



-99 



dx^ 
dr 



99 



.69 



dx^ 



c 9,p,p ko ( d^x 



pI pI 



(21) 

The light surface, which is a generalization of the light cylinder, is obtained by setting ko = 9tt+'2gt^rt+g^^rt'^ 
to be zero (e.g., Znajek 1977; Takahashi et al. 1990). It follows from equation (12) that the acceleration 
electric field is given by E^^ = —{d'^/ds)o,,^,- 

Let us briefly consider equation (19) near the polar cap surface of a nearly aligned rotator. Since 
s « r - r*, 6* < 1, and B^B"^ < B'^, we obtain (Scharlemann, Arons & Fawley 1978, hereafter SAF78; 
Muslimov & Tsygan 1992, hereafter MT92) 



J__9 

J.2 g,^ 



,9£ 
dr 



1 



r2(l — Tg/r) 





2 


'Id/ 




1 d'^-f' 


[d0, 




6 if do if \ 







47r(p- pGj). 



(22) 



Noting that the solid angle element in the metric of magnetic coordinates is given by (to the lowest order in 
0'), 

,B{O,0^,ip^) . 2 



9e.e,del+9^,^J^l = n^iy?^^ {de 



sm' 



we find that the factor 



de, 
de 



* B{s,6*,ip^) 

1 B{s, 6*,, Lp, 



0*dipl) , 



rlBiO,e,,p,) 



(23) 
(24) 



expresses the effect of magnetic field expansion in equation (22). In the same manner, in the general 
equation (19), magnetic field expansion effect is essentially contained in g^'^', 5^**'*, g'^''^', or equivalently, 

in the coefficients of the second-order trans-field derivatives. In what follows, we assume that the azimuthal 
dimension is large compared with the meridional dimension and neglect (p» dependences. 



2.3. Particle Boltzmann Equations 

The second kind equations we have to consider is the Boltzmann equations for particles. At time t, 
position r, and momentum p, the distribution function of positrons (or N- of electrons) obeys the 
following Boltzmann equation, 

^+^.^N±+(^qE+'^xB).^ = S±{t,r,p), (25) 

where v =p/{meT); rrie refers to the rest mass of the electron, q the charge on the particle, and F = 
— {\v\/c)^ the Lorentz factor. In a pulsar magnetosphere, the collision term ^.i. (or 5_) consists of 
the terms representing the appearing and disappearing rates of positrons (or electrons) at r and p per unit 
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time per unit phase- space volume due to pair creation, pair annihilation, and the energy transfer due to IC 
scatterings and synchro-curvature process. 

Imposing a stationary condition 

utilizing V B = 0, and introducing dimensionless particle densities per unit magnetic flux tube such that 
n± = N±/ (p,B /2-Kce), we can reduce the particle Boltzmann equations as 

dn± dpdn± dx dn± 

where the upper and lower signs correspond to the positrons (with charge q = +e) and electrons {q = — e), 
respectively; p = \p\ and 

— =qEi\ COS X (28) 

at c 

dx _ qEii sinx ^ djluB^/^) . 

— = + c^— sm X, (29) 

dt p OS 

ds 

— =ccosx. (30) 

(Introduction of the radiation- reaction force, Psc/c, will be discussed in the next paragraph.) For outward- 
(or inward-) migrating particles, cosx > (or cosx < 0). Since we consider relativistic particles, we 
obtain T = p/{mec). The second term in the right-hand side of equation (29) shows that the particle's 
pitch angle evolves due to the the variation of B (e.g., § 12.6 of Jackson 1962). For example, without 
inward- migrating particles would be reflected by the magnetic mirrors. Using n±, we can express Pe as 

^B f f 

Pe = ^ // bT'+i^^^*^V*^'^^x) - n-{s,9^,ip*,T,x)]drdx. (31) 

The radiation-reaction force due to synchro-curvature radiation is given by (Cheng & Zhang 1996; Zhang & 
Cheng 1997), 

Psc e^r^Q2 , 7 



1 + 37^ ' (32) 



where 



c 12rc V r^Q 



(r-B +Pc)(ccosx/Pc)^ +r-Bw|' 



re V Pc Pc r-B 

r.2 



_rmeC^sinx _ eB 

rB = 5 , wb = (35) 

eB TirieC 

and pc is the curvature radius of the magnetic field line. In the limit of x — * (or pc — > oo), equation (32) 
becomes the expression of pure curvature (or pure synchrotron) emission. 

Let us briefly discuss the inclusion of the radiation-reaction force, Psc/c, in equation (28). Except 
for the vicinity of the star, the magnetic field is much less than the critical value (Bcr = 4.41 x 10^'^ G) 
so that quantum effects can be neglected in synchrotron radiation. Thus, we regard the radiation-reaction 
force, which is continuous, as an external force acting on a particle. Near the star, if r{B/Bcr)smx > 0.1 
holds, the energy loss rate decreases from the classical formula (Erber et al. 1966). If r(B/Bcr)sinx > 1 
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holds very close to the star, the particle motion perpendicular to the field is quantized and the emission is 
described by the transitions between Landau states; thus, equation (28) and (32) breaks down. In this case, 
we artificially put x = 10~^°, which guarantees pure-curvature radiation after the particles have fallen onto 
the ground-state Landau level, avoiding to discuss the detailed quantum effects in the strong-B region. This 
treatment will not affect the main conclusions of this paper, because the gap electrodynamics is governed by 
the pair creation taking place not very close to the star. 

Collision terms are expressed as 

S±{s,d*,ip*,T,x) = - rf£'^7?7c(^7'r,Mc)"±(s,6'*,</'*,r,x) 
J-1 Je^<t 

-I- / d/ic / dri?7fc(ri,r,jUc)n±(s,6i*,</9*,ri,x) 

J-l JTi>T 



+ duc dE. 



^1 



+ ^ ) -B-9±{r,E^,K) 



dT dV B 



(36) 



where refers to the cosine of the collision angle between the particles and the soft photons for inverse- 

Compton scatterings (ICS), between the 7-rays and the soft photons for two-photon pair creation, and 
between the 7-rays and the local magnetic field linos for one-photon pair creation. The function g represents 
the 7-ray distribution function divided by Q.B^/{2-Kce) at energy E-^, momentum k and position r, where 
denotes the polar-cap magnetic field strength. Here, k should be understood to represent the photon 
propagation direction, because and k arc related with the dispersion relation (see next section). Since 
pair annihilation is negligible, we do not include this effect in equation (36). 

If we multiply dV on both sides of equation (36), the first (or the second) term in the right-hand 
side represents the rate of particles disappearing from (or appearing into) the energy interval rUeC^r and 
meC^(r -|- rfr) due to inverse-Compton (IC) scatterings; the third term does the rate of two-photon and 
one-photon pair creation processes. 

The IC redistribution fimction rfl^{E^ ,T , ^i^) represents the probability that a particle with Lorentz 
factor r upscatters photons into energies between E^ and E^ + dE^ per unit time when the collision angle is 
cos~^/Uc. On the other hand, vifQiTi^T^iic) describes the probability that a particle changes Lorentz factor 
from Fj to r in a scattering. Thus, energy conservation gives 

r7?c(r.,r/,Mc) = 7?7c[(r. - Tf)m,c\Vi,iJi,]. (37) 

The quantity ry^c defined by the soft photon fiux dF^/dEg and the Klein-Nishina cross section ctkn as 
follows (HHS03): 

r?7c(i^7>r,/Xc) = (l-/3Mc) 



where (3 = y'l — l/F^ is virtually unity, Q.^ the solid angle of upscattercd photon, the prime denotes the 
quantities in the electron (or positron) rest frame, and E^ = (&i_i -|- bi)/2. In the rest frame of a particle, 
a scattering always takes place well above the resonance energy. Thus, the Klein-Nishina cross section can 
be applied to the present problem. The soft photon flux per unit photon energy E^ [s~^cm~^] is written as 



dFs/dEs and is given by the surface blackbody emission with redshift corrections at each distance from the 
star. 



The differential pair-creation redistribution function is given by 



dr]. 



77 



dr 



poo 

(i;^,r,/ie) = (i-Mc) / 



where the pair-creation threshold energy is defined by 



1 



dFs dap{E^,r,iJ,c) 



dE, 



dT 



I- licE^ 



and the differential cross section is given by 

1-/?, 



dcTp 3 

"dT ~ s'^'^ E. 



CM 



7 



1 + /?cm(2 Mcm) 2/3c]y[(l Mcm)^ 



1 /^cm/^cm (1 /^cmMcm)^ 
ctt refers to the Thomson cross section and the center-of-mass quantities are defined as 

2rmoc2 ~E^ „o . 2{m^c^f 



PcmE^ 



/3, 



CM 



1 



(1 - n^)E,E^ 



(39) 



(40) 



(41) 



(42) 



Since a convenient form of drj^B/dV is not given in the literature, we simply assume that all the 
particles are created at the energy FmeC^ = Ej/2 for magnetic pair creation. This treatment does not affect 
the conclusions in the present paper. 

Let us briefly mention the electric current per magnetic flux tube. With projected velocities, ccosx, 
along the field lines, electric current density in units of Q.B/{2'k) is given by 

jgap(s,6'*) =je{s,0^,) +jion{0*), (43) 

where 

je =jj (n- + n+) cos X dpdx; (44) 

.iion denotes the current density carried by the ions emitted from the stellar surface. Since dp/dt and dx/dt 
in equation (27) depend on momentum variables p and X) Je ^.nd hence jgap does not conserve along the field 
line in an exact sense. 

Nevertheless, jgap is kept virtually constant for s. This is because most of the particles have relativistic 
velocities projected along the magnetic field lines at each point. For example, consider a situation that a 
pair is created inwardly at potion s = si. In this case, the positron will return after migrating a certain 
distance (say, Ss, which is positive). In si — Ss < s < si, the positron does not contribute for the electric 
current, because both the inward and the outward current cancel each other in a stationary situation we 
are dealing with, provided that the projected velocity along the field line is relativistic. Only when the 
positronic trajectory on (,s,pcosx) space becomes asymmetric with respect to the pcosx = axis, owing 
to the synchrotron radiation, which is important if | cosx| <C 1 (see fig. 14 of Hirotani & Shibata 1999a), 
the returning positron has a non- vanishing contribution for the current density at s ~ Si — Ss. In s > Si, 
positronic pitch angle is small enough to give a spatially constant contribution to the current density (per 
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magnetic flux tube). For electrons, it always has an inward relativistic projected velocity and hence gives a 
spatially constant contribution to the current density. In practice, the contribution of the returning particles 
with an asymmetric trajectory around x ~ 90° on the current density, can be neglected when we discuss 
the igap- Thus, we can regard that jgap is virtually conserved even though dp/dt and dx/dt have p and x 
dependences. 



2.4. Gamma-ray Boltzmann Equations 

The third kind equations we have to consider is the Boltzmann equation for 7-rays. In general, the 
distribution function g of the 7-rays with momentum fe obeys the following Boltzmann equation 

^ + c^^-Vg{t,r,k) = S^{t,r,k), (45) 

where = —k^ki; is given by 

S, = - df^, dr ^^P^^'^^' • ff(r, E„ k'/k\ k^/k^) 

/I nOO Q 

duc J dTrilf.{E^,T,Hc)-^n±{s,e^,tp^,T,x) 

P-K POO Q 

+ dx drr]sc{E^,T,x)-s-n±{s,e.,^.,T,x), (46) 
Jo Ji ^* 

where rjsc is the synchro-curvature radiation rate [s~^] into the energy interval between and E^ + dE^ 
by a particle migrating with Lorentz factor F, x the pitch angle of particles. Explicit expression of rysc is 
given by Cheng and Zhang (1996). 

Imposing the stationary condition (26), or equivalently, assuming that g depends on and t bs g = 
g{r, 9,ip — ^It, fe), we obtain 



kf \ dg k-" dg k^ dg 



cttt - ) ^ + CiTTif + cjTjT^ = S^{r, 6, <p, c|fe|. A;'', k^, k^), (47) 



where Cp = ip — fit. To compute A:*, wc have to solve the photon propagation in the curved spacetime. Since 
the wavelength is much shorter than the typical system scales, geometrical optics gives the evolution of 
momentum and position of a photon by the Hamilton-Jacobi equations, 

dkr^ _ _dkt dfc0__9fct 

dX ~ dr' dX ~ 09 ^ ' 

dr_ _ dkt_ d9^ _ dkt_ , , 

dX ~ dkr ' dX^ dke' ^ ^ 

where the parameter A is defined so that cdX represents the distance (i.e., line element) along the ray path. 
The photon energy at infinity kt and the azimuthal wave number — fc^, are conserved along the photon path 
in a stationary and axisymmetric spacetime (e.g., in the spacetime described by eqs. [l]-[3]). Hamiltonian 
kf can be expressed in terms of kr, kg, k^, r, 6 from the dispersion relation fc^fc^ = 0, which is a quadratic 
equation of fc^ (/U = t, r, 0, ip). Thus, we have to solve the set of four ordinary differential equations (48) and 
(49) for the four quantities, k^, k$, r, and 6 along the ray. Initial conditions at the emitting point are given 
by fc*/|fe| = ±B^/\B\, where i = r, 6, (p; the upper (or lower) sign is chosen for the 7-rays emitted by an 
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outward- (or inward-) migrating particle. When a photon is emitted with energy -Eiocai by the particle of 
which angular velocity is ip, it is related with kt and —k^ by the redshift relation, E'locai = {dt/dT){kt + k^ip), 
where dt/dr is solved from {dt / dr)"^ {gtt + 2gt^(fi + gipipifi"^) = 1. To express the energy dependence of g, we 
regard 5 as a function of kt = (i.e., observed photon energy). 

In this paper, in accordance with the two-dimensional analysis of equations (19) and (27), we neglect (p 
dependence of by ignoring the first term in the left-hand side of equation (47). In addition, we neglect the 
aberration of photons and simply assume that the 7-rays do not have angular momenta and put k^ = 0. The 
aberration effects are important when we discuss how the outward-directed 7-rays will be observed. However, 
they can be correctly taken into account only when we compute the propagation of emitted photons in the 
three-dimensional magnetosphere. Moreover, they are not essential when wc investigate the electrodynamics, 
because the pair creation is governed by the specific intensity of inward-directed 7-rays, which are mainly 
emitted in a relatively inner region of the magnetosphere. Thus, it seems reasonable to adopt A;,^ = when 
we investigate the two-dimensional gap electrodynamics. 

We linearly divide the longitudinal distance into 400 grids from s = (i.e., stellar surface) to s = 1.4k7lc, 

and the meridional coordinate into 16 field lines from 0^ = 0"^^^ (i.e., the last-open field line) to 0^ = 0™'" 
(i.e., gap upper boundary), and consider only = plane (i.e., the field lines threading the stellar surface 
on the plane formed by the rotation and magnetic axes). To solve the particle Boltzmann equations (27), we 
adopt the Cubic Interpolated Propagation (CIP) scheme with the fractional step technique to shift the profile 
of the distribution functions n± in the direction of the velocity vector in the two-dimensional momentum 
space (e.g., Nakamura & Yabe 1999). To solve the 7-ray Boltzmann equation (47), on the other hand, we do 
not have to compute the advection of g in the momentum space, because only the spatial derivative terms 
remain after integrating over 7-ray energy bins, which are logarithmically divided from f3i = 0.511 MeV to 
/329 = 28.7 TeV into 29 bins. The 7-ray propagation directions, k^ /k^\ are divided linearly into 180 bins every 
A9y = 2 degrees. Since the specific intensity in ith energy bin at height 9* = 9^ = 6'™''''- (fc/16)(6'™''''-6'™™), 
is given by 

= a/a^ / ' 5[s,^',i^7> {k'/ni] dE^, (50) 
the observed 7-ray energy flux at distance d is calculated as 

where Ay denotes the azimuthal dimension of the gap at longitudinal distance s (= tolc in this paper), 

Azfe the meridional thickness between two field lines with 6^ = 6',^ and ^*jt+i, and i = 1, 2, 3, . . . , 28, 
k — 1, 2, 3, . . . , 15, I — 1,2,3,..., 180. To compute the phase-averaged spectrum, we set the azimuthal width 
of the 7-ray propagation direction, Acp^ to be tt radian. 

Equation (46) describes the 7-ray absorption and creation rate within the gap. However, to compute 
observable fluxes, we also have to consider the synchrotron emission by the secondary, tertiary, and higher- 
generation pairs that are created outside of the gap. If an electron or positron is created with energy FomeC^ 
and pitch angle xo, it radiates the following number of 7-rays (in units of f2-B*/27rce) in energies between 
bi-i and 6j: 

dgi 2-Kce f°° f''^ 1 dW 

-7^ = TTT^ / dt / ——r^dE^, (52) 
dn OB* Jo Jb,_^ E-y dtdE-y ^' ^ ' 

where 

.1^ ^/^^^«--VfM, (53) 



dtdE-y hm^c^ \ Ec 
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"'''^ dt ~ 3 m2c3 ' ^^^> 

F{x) = X / K,/s{Od^; (55) 

is the modified Bessel function of 5/3 order, and Ec = {3h / 4:1t) {eBF'^ sin Xi) / {mec) is the synchrotron 
critical energy at Lorentz factor F. Substituting equations (53) and (54) into (52), we obtain 



dn 0-B* 2/ieSsinxo Ji 



Jbi-i/Ec Jy 



K^/sim- (56) 



Note that we assume that particle pitch angle is fixed at x = xoi because ultra-relativistic particles emit 

radiation mostly in the instantaneous velocity direction, preventing pitch-angle evolution. Once particles 
lose sufficient energies, they preferentially lose perpendicular momentum; nevertheless, such less-energetic 
particles hardly emit synchrotron photons above MeV energies. On these grounds, to incorporate the syn- 
chrotron radiation of higher-generation pairs created outside of the gap, we add {dn/dTf)){dgi/dn)dTo to 
compute the emission of 7-rays in the energy interval [bi_i,bi\ in the right-hand side of equation (45), where 
dn/dVQ denotes the particles created between position s and s + ds in Lorentz factor interval [Tq, Tq + dTo]. 



2.5. Boundary Conditions 

In order to solve the set of partial differential equations (19), (27), and (47) for ^, n±, and g, we must 

impose appropriate boundary conditions. Wc assume that the gap Zower boundary, 9^, = 9^^^, coincides with 
the last open field line, which is defined by the condition that sin 9y/—grrB^ + cos 9^J—geeB^ = is satisfied 
at the light cylinder on the surface t^, = 0. Moreover, we assume that the upper boundary coincides with 
a specific magnetic field line and parameterize this field line with 9^, =- 9f™ . In general, 0™'" is a function 
of however, we consider only (p* = in this paper. Determining the upper boundary from physical 
consideration is a subtle issue, which is beyond the scope of the present paper. Therefore, we treat 6*™'" as 
a free parameter. We measure the trans-field thickness of the gap with 

^max /3inin 

/ixn = * * ■ (57) 

If /im = 1-0, it means that the gap exists along all the open field lines. On the other hand, if /im <C 1, 
the gap becomes transversely thin and 0, derivatives dominate in equation (19). To describe the trans-field 
structure, we introduce the fractional height as 

h = ^. (58) 

Thus, the lower and upper boundaries arc given by h — and h = h^, respectively. 

The inner boundary is assumed to be located at the neutron star surface. For the outer boundary, 
we solve the Poisson equation to a large enough distance, s = lAwhc, which is located outside of the 
light cylinder. This mathematical outer boundary is introduced only for convenience in order that the 
distribution inside of the light cylinder may not be influenced by the artificially chosen outer boundary 
position when we solve the Poisson equation. Since the structure of the outer- most part of the magnetosphere 
is highly unknown, we artificially set E^^ = if the distance from the rotation axis, zu, becomes greater than 
0.90roLc- Under this artificially suppressed £'11 distribution in w > 0.90zuj,Cj we solve the Boltzmann 



-14- 



equations for outward-migrating particles and 7-rays in < s < 1.4ra7Lc- For inward-migrating particles and 
7-rays, we solve only in zj < O.Qojlc- The position of the mathematical outer boundary {IAwlc in this 
case), little affects the results by virtue of the artificial boundary condition, i^y = for > O.QnjLc- On 
the other hand, the artificial outer boundary condition, = for > 0.9ii7lc, affects the calculation of 

outward-directed 7-rays to some degree; nevertheless, it little affects the electrodynamics in the inner part 
of the gap (s < O.Sc^lc), which is governed by the absorption of inward-directed 7-rays. 

First, to solve the elliptic-type equation (19), we impose = on the lower, upper, and inner boundaries. 
At the mathematical outer boundary (s = lAzuhc), we impose d'^ /ds = 0. Generally speaking, the solved 
E\\ = —{d^/ds)s^o under these boundary conditions does not vanish at the stellar surface. Let us consider 
how to cancel this remaining electric field. 

For a super-GJ current density in the sense that Pe — poj < holds at the stellar surface, equation (19) 
gives a positive electric field near the star. In this case, we assume that ions arc emitted from the stellar 
surface so that the additional positive charge in the thin non-relativistic region may bring i^y to zero (for the 
possibility of free ejection of ions due to a low work function, see Jones 1985, Neuhauser et al. 1986, 1987). 
The column density in the non-relativistic region becomes (SAF78) 

where q/m represents the charge-to- mass ratio of the ions and jion the ionic current density in units of 
ilBf/(2Tv). Equating 47rSNR to — (9"I'/<9s),,^o calculated from relativistic positrons, electrons and ions, we 
obtain the ion injection rate jion that c;ancels i?y at the stellar surface. 

For a sub-GJ current density in the sense that pe — PGJ > holds at the stellar surface, \I' increases 
outwards near the star to peak around s = 0.02ii7lc ^ CIOotlC) depending on Ui and pe{s = 0), then 
decrease to become negative in the outer magnetosphere. That is, — (9\I//9s)j5^o < holds in the inner 
region of the gap. In this case, we assume that electrons are emitted from the stellar surface and fill out the 
region where 4" > 0; thus, we artificially put '5 = if \E' > appears. Even though a non-vanishing, positive 
i^ll is remained at the inner boundary, which is located away from the stellar surface, we neglect such details. 
This is because the gap with a sub-GJ current density is found to be inactive and hence less important, as 
will be demonstrated in the next section. 

Secondly, to solve the hyperbolic- type equations (27) and (47), we assume that neither positrons nor 
7-rays are injected across the inner boundary; thus, we impose 

n+(s-,e„r,x) = 0, ff(s'",e,,^^,0^) = O (60) 

for arbitrary 9^, T, < x < 7r/2, E^, and cos{6^. — 6b) > 0, where 9b designates the outward magnetic field 
direction. In the same manner, at the outer boundary, we impose 

n_(s°"',0„r,x)=O, g{s°''\9,,E^,9^)=0 (61) 

for arbitrary 9^, F, 7r/2 < x < tt, Ej, and cos(^-y — ^b) < 0. 

3. Application to the Crab Pulsar 

Since the formulation described in the foregoing section is generic, we specify some of the quantities in 
§ 3.1 before turning to a closer examination in §§3.2-3.8. 
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3.1. Assumptions on Magnetic Field and Soft Photon Field 



First, let us specify the magnetic field. Near the star, we adopt the static (unperturbed by rotation 
and currents) dipole solution obtained in the Schwarzschild space time (e.g., MT92, and references therein). 
That is, in equations (9), we evaluate Fg^p and F^j. as 



/(r)^cose. 



_^d_ r/(r) 

where @ is the angle measured from the magnetic axis, and 

2 / \ 3 

IT* I r \ 

f{r) = -3 



1 r 
2^ 



sin 6, 



(62) 
(63) 

(64) 



At the high altitudes (but within the light cylinder), the open field lines deviate from the static dipole 
to be swept back in the opposite direction of the rotation and bent toward the rotational equator. There are 
two important mechanisms that cause the deviation: Charge flow along the open field lines, and retardation 
of an inclined, rotating dipole. Both of them appear as the first order correction in w/vj-lq expansion 
to the static dipole. To study the former correction, Muslimov and Harding (2005) employed the space- 
charge-limited longitudinal current solved by MT92, and derived an analytic solution of the correction. 
However, if we discard the space-charge-limited-flow of emitted electrons and consider copious pair creation 
in the gap, we have to derive a more general correction formula that is applicable for arbitrary longitudinal 
current distribution. To follow up this general issue further would involve us in other factors than the 
electrodynamics of the accelerator, and would take us beyond the scope of this paper. Thus, we consider 
only the latter correction and adopt the inclined, vacuum magnetic field solution obtained by CRZOO (their 
equations [B2]-[B4]). 

Secondly, we consider how the toroidal current density, J"^ , affects pcj near the light cylinder. In the 
outer magnetosphere, general relativistic effects are negligible; thus, equation (14) becomes 



PGJ = — ^ h 



27rc 



(65) 



Since J'^ is of the order of {CIB/2it){vu/wi,c), the second term appears as a positive correction which is 
proportional to (iu/culc)^ and will become comparable to the first term if ro/tULC ^ 1- Thus, to incorporate 
this special relativistic correction, we adopt 



Pgj = - ■ 



n B 

2irc 



1 



(66) 



where the constant k is of the order of unity. For example, CHR86a adopted k = \. Even though a larger 
value of K is preferable to reproduce a harder curvature spectrum above 5 GeV and a larger secondary 
synchrotron flux around 100 McV, wc adopt a conservative value k ~ 0.5 in the present paper. 

Thirdly, we have to specify the differential soft photon flux, dFg/dEs, which appears in equations (38) 
and (39). As the possible soft photon fields illuminating the gap, we can consider the following three com- 
ponents in general: 
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(1) Photospheric emission from the hot surface of a cooHng neutron star. For simphcity, we approximate 
this component with a black body spectrum with a single temperature, kTg. We assume that this component 
is uniformly emitted from the whole neutron star surface. 

(2) Thermal soft X-ray emission from the neutron star's polar cap heated by the bombardment of relativistic 
particles streaming towards the star from the magnctosphcrc. Since we consider a young pulsar in this paper, 
this component is negligible compared to the first component. 

(3) Non-thermal, power-law emission from charged relativistic particles created outside of the gap in the 
magnetosphere. The emitted radiation can be observed from optical to 7-ray band. Since the non-thermal 
emission will be beamed away from the gap, we assume that this component does not illuminate the gap. The 
major conclusions in this paper will not be affected by this assumption, except that the pair creation would 
increase to suppress the potential drop and hence the 7-ray luminosity if this component illuminates the gap. 

We apply the scheme to the Crab pulsar, adopting four free parameters, cci, /x, kTg, and /im- Other 

quantities such as gap geometry on the poloidal plane, exerted and potential drop, particle density and 
energy distribution, as well as the 7-ray flux and spectrum, are uniquely determined if we specify these four 
parameters. 

In the next section, we consider transversely thin and thick cases in § 3.2 and §§ 3.3-3.8, respectively. 

3.2. Sub-GJ current solution: Traditional Outer-gap Model 

To begin with, let us consider the magnetic inclination ai = 70°, which is more or less close to the value 
(65°) suggested by a three-dimensional analysis in the traditional outer gap model (CRZOO). We adopt 
kTs — 100 eV as the surface blackbody temperature, which is consistent with the observational upper limit, 
fcTs < 180 eV (Tennant et a. 2001). In §§ 3.2-3.4, we adopt /u = 4.0 x lO^^G cm^ as the magnetic dipole 
moment, which gives B^, — 1.46 x lO^'^G (cq. [10]), assuming = 10® cm and M = 1.4M0. If we evaluate /i 
from the spin-down luminosity E = 2r2''/x^/3c^, we obtain /x = 3.8 x 10"^" G cm"* for E = 4.46 x 10^* ergs s~^ 
(e.g., Becker & Triimper 1997). The dependence of the solution on kTs, fi, and will be discussed in 
§§ 3.4-3.6. 

Examine a sub-GJ current solution, which is defined by je < \pGj / {^B /2ttc)\s=o ~ (1 — u}/fl)Bz*/Bt,, 
where B^* refers to the surface magnetic field component projected along the rotation axis: the right-hand 
side is evaluated at s = and ^* = 6'^'^^. To this aim, we consider a transversely thin gap, hm = 0.047. The 
solution becomes similar to the vacuum one obtained in the traditional outer-gap model model (CHR86a), 
as the left panel of figure 3 indicates. In this figure, we present E\\{s,h) at discrete height h ranging from 
2/ini/16; 5/ini/16; 8ft-m/16, llft-m/lB, 14/ini/16, with dashed, dotted, solid, dash-dot-dot-dot, and dash-dotted 
curves, respectively; they are depicted in the right panel with a larger hm (= 0.200) for illustration purpose. 
For one thing, for such a transversely thin, nearly vacuum gap, the inner boundary is located slightly inside 
of the null surface. What is more, maximizes at the central height, h = /im/2, and remains roughly 
constant in the entire region of the gap. The solved E^^ distributes almost symmetrically with respect to the 
central height; for example, the dashed and dash-dotted curves nearly overlap each other. Similar solutions 
are obtained for a thiner gap, /im < 0.047. 
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3.3. Super-GJ current solution: Hybrid gap structure 

Next, let us consider a thicker gap, /im = 0.048. In this case, je becomes comparable or greater than 
\pG3 /{^B /2ttc)\s=o in the upper half region h > of the gap, deviating the solution from the vacuum 

one. In the left panel of figure 4, we present E^^{s,h) at five discrete height h in the same way as figure 3. 
It follows that is screened by the discharge of created pairs in the inner-most region (s < O.Szui^c) in 
the higher latitudes {h > hra/2). For example, E\\ at h = T/im/S (dash-dotted curve) deviates from the 
unscreened solution at the lower latitudes h = /im/8 (dashed). On the other hand, along the lower field lines 
{h < hni/2), i?||(s) is roughly constant as in the traditional outer-gap model. 

Let us further consider a thicker gap, hm = 0.060. The right panel of figure 4 shows that E\^{s,h) is 
substantially screened than the marginally super-GJ case, hm = 0.048 (left panel). Because the gap transfield 
thickness virtually shrinks in s < 0.5whC due to screening in the higher latitudes, S|| in the lower latitudes 
also decreases compared to smaller hm cases, as the dashed lines in the left and right panels indicate. 

To understand the screening mechanism, it is helpful to examine the Poisson equation (22), which is a 
good approximation in s < OAwj^c- In the transversely thin limit, it becomes 

1 l^W'''*~4,(.-«,)»?5Mf^-^), (67) 



(l-rg/r)r2 dO^ ''^^•^ c \nB/2nc nB/2nc 

Since we are interested in the second-order 6^ derivatives, this equation is valid not only for a nearly aligned 
rotator, but also for an oblique rotator. We could directly check it from the general equation (13). Factoring 
out the magnetic field expansion factor, B{r)/B^, from the both sides, we obtain 

* » (' - t) ^ (nlT^ - Tm^^ -^(«- - »"'"'<""- - «•)■ 

We thus find that E\\ = —d^/ds is approximately proportional to —d{p/B — pQj/B)/ds. It is, therefore, 
important to examine the two-dimensional distribution of p/B and pGj/B to understand E\\{s,h) behavior. 

In figure 5, we present p/(QS/27rc), pe/{flB/2iTc), and pgj/{^B/2-!tc), as the solid, dash-dotted, and 
dashed curves, at nine discrete magnetic latitudes, ranging from h = (4/16)/im, (5/16)/im, ■ . ■, to (12/16)/im 
for the same parameters as the right panel of figure 4. If there is a cold-field ion emission from the star, the 
total charge density (solid curve) deviates from the created charge density (dash-dotted one). It follows that 
the current is sub-GJ for h < (5/16)/im = 0.0187 and super-GJ for h > (6/16)/iin = 0.0225. Along the field 
lines with super-GJ current, p^ — pqj becomes negative close to the star. This inevitably leads to a positive 
i^ll , which extracts ions from the stellar surface. In this paper, we assume that the extracted ions consist of 
protons; nevertheless, the conclusions are little afii'ected by the composition of the extracted ions. 

In the outer region, p/B levels off in s > O.Scclc for h > /i.m/2. Since pcj/B becomes approximately 
a linear function of s, E^^ remains nearly constant in s > O.Sotlc, in the same manner as in traditional 
outer-gap model. In the inner region, on the other hand, inward-directed 7-rays propagate into the convex 
side due to the field line curvature, increasing particle density exponentially with h. As a result, the lower 
part (i.e., smaller h region) becomes nearly vacuum. For example, at h = /im/4 = 0.015 (top left panel), 
positive pe — pgj leads to a negative S|| at the stellar surface, inducing no ion emission. Even though the 
created current is sub-GJ at h = 0.0187, ions are extracted from the surface. This is because the negative 
Pes = P — PGJ in the higher latitudes h > 0.0225 cancels the relatively small positive pes along h = 0.0187 
to induce a positive S|| at the stellar surface. Such a two-dimensional effect in the Poisson equation is also 
important in the higher altitudes (O.Iiulc < s < O.Sojlc) along the higher-latitude field lines(/i > 0.0225). 
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Outside of the null surface, s > O.OQnjLc, there is a negative pcs in the sub-GJ current region {h < 0.015). 
This negative Peff works to prevent E^^ from vanishing in the higher latitudes, where pair creation is copious. 
However, the created pairs discharge until vanishes, resulting in a larger gradient of p than that of pcj in 
the intermediate latitudes in 0.0225 < h < 0.0262. In the upper half region (0.03 < h < hm = 0.06), dp/ds 
does not have to be greater than 9/>gj /ds, in order to screen E\\^ . 

In short, the gap has a hybrid structure: The lower latitudes (with small h) arc nearly vacuum having 
sub-GJ current densities and the inner boundary is located slightly inside of the null surface, because > 
region will be filled with the electrons emitted from the stellar surface. The higher latitudes, on the other 
hand, are non-vacuum having super-GJ current densities, and the inner boundary is located at the stellar 
surface, extracting ions at the rate such that their non-relativistic column density at the stellar surface 
cancels the strong induced by the negative p — pa of relativistic electrons, positrons, and ions. The 
created pairs discharge such that E\^ virtually vanishes in the region where pair creation is copious. Thus, 
in the intermediate latitudes between the sub-GJ and super-GJ regions, dp/ds > poj/ds holds. 

Even though the inner-most region of the gap is inactive, general relativistic effect (space-time dragging 
effect, in this case) is important to determine the ion emission rate from the stellar surface. For example, at 
h = /i.m/2 for = 0.600 (i.e., the central panel in fig. 5), jion is 69% greater than what would be obtained 
in the Newtonian limit, pcj = ■ B/2-kc. This is because the reduced |/9gj| near the star (about 15% 
less than the Newtonian value) enhances the positive E\\ , which has to be canceled by a greater ion emission 
(compared to the Newtonian value). The current, jion, is adjusted so that Ipcff | may balance with the trans- 
field derivative of \E' near the star. The resultant jpeffl becomes small compared to \pGj\, in the same manner 
as in traditional polar-cap models, which has a negative E\\ with electron emission from the star. Although 
the non-relativistic ions have a large positive charge density very close to the star (within 10 cm from the 
surface), it cannot be resolved in figure 5. Note that the present calculation is performed from the stellar 
surface to the outer magnetosphere and does not contain a region with E\\ < 0. It follows that an accelerator 
having E^^ < (e.g., a polar-cap or a polar-slot-gap accelerator) cannot exist along the magnetic field lines 
that have an super-GJ current density created by the mechanism described in the present paper. 

It is worth examining how E^^ changes with varying h^. In figure 6, we present E\\{s,h) at the central 
height h — /im/2. In the left panel, the dotted, solid, dashed, and dash-dotted curves correspond to = 
0.047, 0.048, 0.060, and 0.100, respectively, while in the right panel, the dash-dotted, dash-dot-dot-dotted, 
solid, and dashed ones to 0.100, 0.160, 0.200, and 0.240, respectively. It follows that the inner part of 
the gap becomes substantially screened by the discharge of created pairs as ft-m increases. It also follows 
that the maximum of iSy increases with increasing hm for hm < 0.2, because the two-dimensional screening 
effect due to the zero-potential walls becomes less important for a larger hm- To solve particle and 7-ray 
Boltzmann equations, we artificially put S|| = in ti7 > O.QrajLC, or equivalently in s > l.lrajLC for on = 70°, 
as mentioned in § 2.5. 

The created current density, je, is presented in figure 7, as a function of h. The thin dashed line 
represents \pQ_-f / {ilB /2'!tc)\s=o; if jc appears below (or above) this line, the created current is sub- (or super- 
) GJ along the field line specified by h. The open and filled circles denote the lowest and highest latitudes 
that are used in the computation. For = 0.047, the solution (dotted curve) is sub-GJ along all the 
field lines; thus, screening due to the discharge is negligible as the dotted curve in the left panel of figure 6 
shows. As hm increases, the solution becomes super-GJ from the higher latitudes, as indicated by the solid 
(V = 0.048), dashed (/im = 0.060), dash-dotted {h^ = 0.100), and dash-dot-dot-dotted (V = 0.160) curves 
in figure 7. As a result, screening becomes significant as hm increase, as figure 6 shows. This screening of E'n 
has a negative feed back efii'ect in the sense that je is regulated below unity. Even though it is not resolved 
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in figure 6, in the lower latitudes, jc grows across the gap height exponentially, as CHR86a suggested. For 
example, je = 3.0 x 10'^°, 2.1 x 10-^ 2.1 x 10"^, 1.0 x 10"^ 4.3 x 10"^, 1.0 x lO'^, 4.6 x 10-^ and 
1.3 x 10~^ at h/hm = 1/16, 2/16, 3/16, . . ., 8/16, respectively, for hm = 0.048 (solid curve). This is because 
the pair creation rate at height h is proportional to the number of 7-rays that are emitted by charges on all 
field lines below h. 

Let us now turn to the emitted 7-rays. Figure 8 shows the phase-averaged 7-ray spectrum calculated 
for three different h^^s. Open circles denote the pulsed fluxes detected by COMPTEL (below 30 MeV; 
Ulmer ct al. 1995, Kuiper et al. 2001) and EGRET (above 30 MeV; Nolan et al. 1993, Fierro et al. 1998), 
while the open square does the upper limit obtained by CELESTE (de Naurois et al. 2002). It follows that 
the sub-GJ (i.e., traditional outer-gap) solution with = 0.047 predicts too small 7-ray flux compared 
with the observations. (Note that in traditional outer-gap models, particle number density is assumed to 
be the Goldreich-Julian value, while £^|| is given by the vacuum solution of the Poisson equation, which 
is inconsistent.) The maximum flux, which appears around GeV energies, does not become greater than 
lO^'^ Jy Hz for any sub-GJ solutions, whatever values of ofj, /i, and fcTg we may choose. Thus, we can rule 
out the possibility of a sub-GJ solution for the Crab pulsar. 

As /im increases, the increased i^y results in a harder curvature spectrum, as the solutions corresponding 
to /im = 0.100 and 0.200 indicate. As will be discussed in § 4.1, the problem of insufficient 7-ray fluxes may 
be solved if we consider a three-dimensional gap structure. However, the secondary synchrotron flux emitted 
outside of the gap is too small to explain the flat spectral shape below 100 MeV. The 7-ray spectrum is nearly 
unchanged for 0.2 < /im < 0.3. For hm > 0.3, the 7-ray flux tends to decrease, because the -E||(s) peaks 
outside of the artificial outer boundary, rsinO = 0.9k7lC) which corresponds to s = l.lnJLC for ai = 70°. For 
/),,„ > 0.4, the gap virtually vanishes because of the discharge of copiously created pairs; in another word, 
the gap is located outside of rsin6' ~ 0.9i37lc- On these grounds, we cannot reproduce the observed flat 
spectral shape if we consider kTg = 100 eV, /x = 4.0 x 10"^° G cm"^, and ai = 70°, no matter what value of 
hm is adopted. Therefore, in the next three subsections, we examine how the solution changes if we adopt 
different values of kTg, fj,, and ai. 

3.4. Dependence on Surface Temperature 

In the same manner as in § 3.3, we calculate E\\{s,h = for kT^ = 150 eV to find that their 

distribution is similar to fcTg = 100 eV case (i.e., fig. 6). For example, sub/super-GJ current solutions 
are discriminated by the condition whether is greater than 0.048 or not, and the maximum of i^y is 
7.2 X 10^ V m~^ for 0.20 < < 0.24. Other quantities such as the particle and 7-ray distribution fimctions 
are also similar. Therefore, we can conclude that the solution is little subject to change for the variation of 
kTg, even though the photon-photon pair production rate increases with increasing fcTg. This is due to the 
negative feedback effect, which will be discussed in § 4.2. 

3.5. Dependence on Magnetic Moment 

Let us examine how the solution depends on the magnetic moment, /U. In figure 9, we present £^||(s, h = 
/im/2) for seven discrete h^s with /x = 6.0 x 10^° G cm^; in the left panel, the dotted and solid curves 
correspond to h = 0.039 and 0.041, respectively, instead oi h = 0.047 and 0.048 in figure 6. It follows that 
the exerted iSy is greater than the case of /x = 4.0 x 10^° G cm"^ (fig. 6), because pcj increases 1.5 times. The 
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distance along field line 



Fig. 3. — Left: The field-aligned electric field of a sub-GJ current solution at five discrete heights (see right 
panel) for ai = 70° and /im = 0.047. The abscissa indicates the distance along the field line from the star 
in the unit of the light-cylinder radius. The null surface position at the height h = /im/2 is indicated by the 
down arrow. Right: Magnetic field lines on the poloidal plane in which both the rotational and magnetic 
axes reside. Instead of /im = 0.047, ~ 0.200 is adopted for clarity. The thick solid curves denote the 
lower and upper boundaries, while the thin dashed, dotted, solid, dash-dot-dot-dotted, and dash-dotted ones 
give the same h/h^ values as those in the left panel; they are /i//im = 2/16, 5/16, 8/16, 11/16, and 14/16, 
respectively. 
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Fig. 5. — Total (solid), created (dash-dotted), and Goldreich- Julian (dashed) charge densities in 
^lB{s, h) / {271 c) unit, for ai — 70° and /im = 0.060 at nine transfield heights, h. If there is an ion emis- 
sion from the stellar surface, the total charge density deviates from the created one. 
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distance along field line 



Fig. 6. — Field-aligned electric field at /i = as a function of s/tzjlc for ^JL = 4.0 x 10"^" G cm'^, 

kTs — 100 eV, and a-i — 70°. Left: The dotted, solid, dashed, and dash-dotted curves corresponds to 
h = 0.047, 0.048, 0.060, and 0.100; right The dash-dotted, dash-dot-dot-dotted, solid, and dashed curves 
corresponds to h = 0.100, 0.160, 0.200, and 0.240. 




fractional tronsfield height 

Fig. 7. — Created current density (in unit of ilB/2TT) as a function of the transfield thickness h for 
^ = 4.0 X 10^" G cm^, kTs = 100 eV, and = 70°. The dotted, sohd, dashed, dash-dotted, and dash-dot- 
dot-dotted curves corresponds to h — 0.047, 0.048, 0.060, 0.100, and 0.160, respectively. 
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negative feedback effect cannot cancel the increase of ^, unlike the increase of kTg, because the right-hand 
side of equation (19) is more directly affected by the variation of through pcj than by the variation of kTg 
through pair creation. 

We also present the 7-ray spectrum for three different /x's in figure 10. It follows that both the peak 
energy and the flux of curvature 7-rays (around GeV energies) increase with increasing /x. This is because 
PgJ) and hence increases with /x. It also follows that the secondary synchrotron flux (below 100 MeV) 
increases with fj,, because the magnetic field strength increases in the magnetosphere. We find that a larger 
magnetic dipole moment, /x > 6 x 10"^° Gem'' is preferable to explain the observed pulsed flux from the Crab 
pulsar. 

If we adopt /x = 8 x 10^° G cm^, which is about twice larger than the dipole deduced value, the spectral 
shape becomes more consistent compared with smaller /i cases. We should notice here that the moment of 
inertia, /, have to be large in this case. For example, if we assume a pure magnetic dipole radiation, the 
spin-down luminosity becomes Lsd = {2sin^ ai/3){iJ.'^n^ /c^) — 1.7 x lO''^ ergs s^^ for ai = 70°. Equating 
—E = —Iflfl with this Z/sd, we obtain / = 3.9 x 10^^ g cm^, which is consistent with the limit (/ > 
3.04 X 10^^ g cm^) derived from the consideration of energetics of the Crab nebula (Bcigcr & Hacnsel 2002). 
However, solving the time-dependent equations of force-free electrodynamics, Spitkovsky (2006) derived 
-^SD ~ (1-l-sin^ ai){iJ.'^^l^/c^), which gives w 5.6 x 10^® ergs for ai = 70°. This large spin-down luminosity 
results in 7 Ri 1.2 X 10^^ g cm^, which is too large even compared with those obtained for stiff equation 
of state (Serot 1979a, b; Pandharipandc & Smith 1975). Thus, it implies either that /i should be less 
than 8 x lO'^" Gem'' or that the deduced magnetospheric current is too large when the relationship Lsd ~ 
(1 + sin^ ai)(/x^O^/c'^) is derived (for a discussion of the magnetospheric current determination, see section 
4 of Hirotani 2006). 

3.6. Dependence on Magnetic Inclination 

Let us examine how the solution depends on the magnetic inclination angle, ai. Figure 11 shows the 
7-ray spectra for /z = 6 x 10^° G cm^, kTg = 100 eV, = 0.200, with three different inclination, ai = 50°, 
60°, and 80°, in the same manner as in figure 8. The flux is averaged over the meridional emission angles 
between 44° and 58° (sohd), 58° and 72° (dashed), 72° and 86° (dash-dotted), 86° and 100° (dotted), 
100° and 114° (dash-dot-dot-dotted), from the magnetic axis on the plane in which both the rotational and 
magnetic axes reside. We find that the 7-ray flux reaches a peak of ~ 2 x 10^^ Jy Hz around 2 GeV and 
that this peak does not strongly depend on a^. This is because; the pair creation efficiency, which governs 
the gap electrodynamics, crucially depends on the distance from the star, which has a small dependence on 
ai. However, we also find that the flux tends to be emitted into larger meridional angles from the magnetic 
axis (i.e., from outer regions) for smaller ai. The reasons are fourfold: 

(1) pgj decreases with decreasing ai at a fixed s; as a result, the null surface appears at larger s for smaller 
ai. 

(2) pGj(s) = — PGj(O) is realized at larger s for smaller Ui. 

(3) In the region where pcj < — Pgj(O) holds, E^^ is substantially screened by the discharge of created pairs. 
(Compare fig. 5 and the right panel of fig. 4.) 

(4) The unscreened E\\ tends to appear at larger s for smaller ai, resulting in a 7-ray emission which 
concentrate in larger meridional angles. 

We can alternatively interpret the explanation above as follows: 
(a) |/9gj(0)| increases with decreasing ai. 



- 24 - 



distance— 2 kpc distance— 2 l<pc distance— 2 l<pc 
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Fig. 8. — Calculated phase-averaged spectra of the pulsed, outward-directed 7-rays for ai — 70°, kTg — 
100 eV, and fi — 4.0 x 10"^° G cni"^, with three different gap thickness, /im- The flux is averaged over the 
meridional emission angles between 44° and 58° (solid), 58° and 72° (dashed), 72° and 86° (dash-dotted), 
86° and 100° (dotted), 100° and 114° (dash-dot-dot-dotted), from the magnetic axis on the plane in which 
both the rotational and magnetic axes reside. 
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distance along field line distance along field line 

Fig. 9. — Same figure as figure 6 but with /i = 6.0 x 10^" G cm'^; in the left panel, the dotted and solid curves 
correspond to h = 0.039 and 0.041, respectively. For other two curves in the left panel and four curves in 
the right panel, h takes the same values as in figure 6. 



- 25 - 



(b) The created current density, jo, which is greater than c|/Ogj(0)| for a super-GJ solution, increases with 
decreasing aj. For example, we obtain jo = 0.62, 0.67, 0.74, and 0.77 for ai = 80°, 70°, 60°, and 50°, 
respectively, at /i = = 0.100 with ^ = 6 x 10^° G cm^, and kT^ = 100 cV. 

(c) A larger jc results in a larger injection of the discharged positrons and the emitted ions into the strong- 
region from the stellar side. 

(d) A larger injection of positive charges from the stellar side shifts the gap outwards by the mechanism 
discussed in § 2 of Hirotani and Shibata (2001a). 

For a smaller inclination, aj < 40°, the gap is located in tn > 0.9tULc; that is, no super-GJ solution exists. 
Therefore, the observed 7-ray flux cannot be explained by the present theory, if the magnetic inclination is 
constrained to be less than 40° by some other methods. 



3.7. Particle Distribution Functions 

Let us examine how the particle distribution function evolves at different positions. Figure 12 represents 
the evolution of positronic distribution function from s — 0.85n7Lc (dashed curve), CQOtolc (dash-dotted), 
0.95wlc (dash-dot-dot-dotted), l.OOtULc (solid), to LSOzhlc (dotted). It shows that the positrons are 
injected into the strong- region (s > O.SStULc) with energies 10^ < F < 3 x 10^ because of the acceleration 
by the small-amplitude, residual in s < O.BStnLc (solid curve in the right panel of fig. 9). The positrons are 
accelerated outwards to attain F ^ 4 x 10^ at s ~ l.OtULc- They are subsequently decelerated by curvature 
cooling in nj > O.Owlc (or equivalently s > I.Iccjlc for a-i — 70°), where we artificially put = 0, to 
escape outwards with F ^ 10^ at s ~ I.Siulc- There is a small population of the positrons that have 
upscattered surface X-rays to possess smaller energies than the curvature-limited positrons. For example, 
within the gap (s < LIcclc), the dash-dotted, dash-dot-dot-dotted, and solid curves have the broad, lower- 
energy component, which connects with the curvature-limited peak component. However, in s > I.Iwlc, 
we artificially put = 0; as a result, the positrons that have lost energies by ICS cannot be re-accelerated, 
forming a separate component in F < 4 x 10^ from the curvature-limited peak component, as the dotted 
curve shows. The upscattered photons obtain several TeV energies; however, they are totally absorbed by 
the strong magnetospheric infrared radiation field. Therefore, we depict only the photon energies below 
100 GeV in figures 8, 10, and 11. 
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Fig. 10. — Same figure as figure 8 for ai = 70°, kTs = 100 eV, and hm — 0.200, with diff"ereiit dipole moment, 
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Since most of the pairs are created inwards, positrons return to migrate outwards by the small-ampHtude 
in s < 0.85ci7lc (for fj, = 6.0 x 10"^" G cm'^ and ai = 70°), losing significant transverse momenta via 
synchrotron process to faU onto the ground-state Landau level in strong B region. Thus, their emission in 
strong- region is given by a pure-curvature formula. 

Next, we consider the distribution function of electrons. In figure 13 we present their evolution along 
the field line h = /im/2 — 0.100 from s — O.OOtulc (dash-dot-dot-dotted curve), 0.40tULc (dotted), 0.20tULc 
(dash-dotted), O.OSgjlc (dashed), to (solid). Left panel shows the energy spectrum. Electrons created 
in s < 0.85n7Lc cannot be accelerated by £^|[ efficiently; thus, their energy spectrum becomes broad as 
the dotted, dashed-dotted, and dashed curves indicate (i.e., particle Lorentz factors do not concentrate at 
the curvature-limited terminal value). From s = O.OStULc to s = 0, electrons are decelerated via synchro- 
curvature radiation and non-resonant inverse-Compton scatterings. Finally, they hit the stellar surface with 
F < 3 X 10^. Assuming that the azimuthal gap width is tt radian, we obtain Lpc = 2.6 x lO'^^ ergs s~^ as the 
heated polar-cap luminosity. Thus, the X-ray emission due to the bombardment is negligible, compared with 
the total soft-X-ray luminosity (0.1-2.4keV) of 7.6 x lO^"* ergs s-i(Arix/ster) (e.g., Becker & Triimper 1997), 
where Afix refers to the emission solid angle. 

Let us see the pitch angle evolution of the electrons, which is presented by the right panel of figure 13. 
In the outer part of the gap, the electrons are created by the collisions between the outward-directed 7-rays 
and the surface X-rays. Thus, created electrons have outward momenta initially, then return by the positive 

losing their perpendicular momentum substantially via synchrotron radiation. Thus, at s = 0.90ci7LC; 
the dash-dot-dot-dotted curve show that their pitch angles, are less than 8 x 10~^. However, most of the 
pairs are created by the inward-directed 7-rays in s < 0.6wlc', thus, electrons have initial inward momenta to 
migrate inwards by the small-amplitude, residual . Since such inward-created electrons do not change their 
migration direction, their pitch angles are greater than those of the outward-created ones, as the dotted curve 
demonstrates. As the electrons migrate inwards, they lose perpendicular momenta via synchro-curvature 
radiation in the strong B field, reducing their pitch angles, as the dotted, dash-dotted, and dashed curves 
indicate. 

We should point out that a pure curvature formula cannot be applied to the electrons. For example, 
at s = OAwi^c, the dotted curve in figure 13 demonstrates that electrons have 10'* < F < 10^ and 10~* < 
sinx < 10~^'^. Noting that we have B = 3.0 x 10^ G at s = 0.4-c<7lc, we find that newly created electrons, 
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Fig. 11.— Same figure as figure 8 for kTs = 100 eV, fi = 6.0 x 10^° G cm^, and /im = 0.200 with different 
magnetic inclination angle, ai. 
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which have lower-energies (F < 10^) and larger pitch angles (sin% > 10~^), emit synchro-curvature radiation, 
rather than pure-curvature one, as figure 14 shows. At s = 0.2wlc, we have B = 3.6 x 10^ G; thus, the dash- 
dotted curve shows that the pure-curvature formula is totally inapplicable as we consider a smaller distance 
from the star. The only exception is the inner-most region (s < 5r* = O.OSlroLC and B > 5 x 10^° G), 
where electrons suffer substantial dc-cxcitation via synchrotron radiation, falling at last onto the ground- 
state Landau level. In this region, electrons emit via pure curvature radiation. We thus artificially assume 
sinx = 10~^°, which guarantees pure curvature emission, and do not depict the solid curve in the right 
panel of figure 13. Since pair creation and the resultant screening of E^^ is governed by the inward-directed 
7-ray flux and spectrum, it is essential to adopt the correct radiation formula by computing the pitch-angle 
evolution of inward-migrating particles. 



3.8. Formation of Magnetically Dominated Wind 

Let us finally consider the magnetic dominance within the light cylinder. First, introduce the magneti- 
zation parameter, 

a = ^ (Pmec^ /// (iV+ + N_)d^p + Fionmionc" /// VViond^p 



4n 



e(B/2)n7Lc 



(69) 



lon^ ion"''ion»- 



where iVion, Lion, refers to the distribution function and the averaged Lorentz factor of the ions. Evaluating 
equation (69) at the light cylinder, and noting Fion^TiionC^ = leA^fgapl, where |eA^gap| ~ (B/4)/i^a7LC refers 
to the potential drop in the gap (see eq. [68]), we obtain 

1 2 

"^"^ ~ /4 ' iermeCV|eA*gap|+iion" ^^^^ 

Recalling jion < je < 1, FmoC^ < |eA4'gap|, and hm < 1, we can conclude that the magnetic energy fiux 
is always greater than the particle kinetic energy flux at the light cylinder, regardless of the species of the 
accelerated particles, of the sign of i^n, or of the gap position (i.e., whether inner or outer magnetosphere) . 
Only when the gap is formed along most of the open field lines (i.e., /im ~ 1), cl can be of the order of unity. 
Since the second factor in equation (70) does not change significantly for different values of parameters 
except for h^, ai^ solely depends on /im. Substituting jion ^ 0.3, jc ^ 0.6, FmcC^/jeA^'gapl ^ 0.25 in 
(3/16)/ini <h< (13/16)/ini for a; = 70°, /i = 6 x 10'™ G cm-\ and h^-a = 0.20, we obtain (Tl ^ 110. 

In short, along the open field lines threading the gap, Poynting flux dominates particle kinetic energy 
flux by a factor of 100 at the light cylinder. Positrons escape from the gap with F lO'', while ions with 
V ~ 1 0^ 

ion ^ J-W . 



4. Summary and Discussion 

To conclude, we investigated the self-consistent electrodynamic structure of a particle accelerator in the 
Crab pulsar magnetosphere on the two-dimensional surface that contains the magnetic field lines threading 
the stellar surface on the plane in which both the rotation and magnetic axes reside. We regard the trans- 
field thickness, /im, of the gap as a free parameter, instead of trying to constrain it. For a small /im, the 
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Lorentz factors of positrons 

Fig. 12. — Energy spectrum of positrons &t s — 0.85ti7Lc (dashed), 0.90roLc (dash-dotted), 0.95n7Lc (dash- 
dot-dot-dotted), l.OOtULc (solid), and I.SOwlc (dotted), for a; = 70°, ^ = 6.0 x 10^° G cm^ kT^ = 100 eV, 
and /im = 0.20. 




Lorentz factors of electrons sin(pitch angles) 



Fig. 13. — Distribution function of electrons at s = 0.9tt7LC (dash-dot-dot-dotted), OAvJhc (dotted), 0.2zuhC 
(dash-dotted), O.OSwlc (dashed), and (solid), for ai = 70°, ^ = 6.0 x lO^o G cnl^ kT^ = 100 eV, and 
/i,n = 0.20. Left panel shows the Lorentz factor dependence, while the right one shows the pitch angle 
dependence. 
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created current density, je, becomes sub-Goldreich- Julian, giving the traditional outer-gap solution but with 
negligible 7-ray flux. However, as hm increases, je increases to become super-GJ, giving a new gap solution 
with substantially screened acceleration electric field, S|| , in the inner part. In this case, the gap extends 
towards the neutron star with a small-amplitude positive which extracts ions from the stellar surface. 
It is essential to examine the pitch-angle evolution of the created particles, because the inward-migrating 
particles emit 7-rays, which governs the gap electrodynamics through pair creation, via synchro-curvature 
process rather than pure-curvature one. The resultant spectral shape of the outward-directed 7-rays is 
consistent with the existing observations; however, their predicted fluxes appear insufficient. The pulsar 
wind at the light cylinder is magnetically dominated: Along the field lines threading the gap, magnetization 
parameter, ul, is about 10^. 

4.1. How to obtain sufficient gamma-ray flux 

The obtained 7-ray fluxes in the present work, are all below the observed values. Nevertheless, this 

problem may be solved if we extend the current analysis into a three-dimensional configuration space. As 
mentioned in § 2.4, we assume kip = and neglect the aberration of photon emission directions. However, in 
a realistic three-dimensional pulsar magnetosphcre, 7-rays will have angular momenta and may be emitted 
in a Hmited sohd angle as suggested by the caustic model (e.g., Dyks & Rudak 2003; Dyks, Harding, & 
Rudak 2004), which incorporate the effect of aberration of photons and that of time-of- flight delays. In 
particular, in the trailing peak of a highly inclined rotator, photons emitted at different altitudes s will 
be beamed in a narrow solid angle to be piled up at the same phase of a pulse (Morini 1983; Romani & 
Yadigaroglu 1995), resulting in a 7-ray flux which is an order of magnitude greater than the present values. 
Thus, the insufficient 7-ray fluxes does not suggest the inapplicability of the present method. It is noteworthy 
that the meridional propagation angles of the emitted photons (e.g., different curves in figs. 8, 10, and 11) 
can be readily translated into the emissivity distribution in the gap as a function of s. Therefore, in this 
work, we do not sum up the 7-ray fluxes emitted into different meridional angles taking account of the 
aberration of light. 

4.2. Stability of the Gap 

Let us discuss the electrodynamic stability of the gap, by considering whether an initial perturbation 
of some quantity tends to be canceled or not. In the present paper, we consider that the soft photon field 
is given and unchanged when gap quantities vary. Thus, let us first consider the case when the soft photon 
field is fixed. Imagine that the created pairs are decreased as an initial perturbation. It leads to an increase 
of the potential drop due to less efficient screening by the discharged pairs, and hence to an increase of 
particle energies. Then the particles emit synchro-curvature radiation efficiently, resulting in an increase of 
the created pairs, which tends to cancel the initial decrease of created pairs. 

Let us next consider the case when the soft photon field also changes. Imagine again that the created 
pairs are decreased as an initial perturbation. It leads to an increase of particle energies in the same 
manner as discussed just above. The increased particle energies increase not only the number and density of 
synchro-curvature 7-rays, but also the surface blackbody emission from heated polar caps and the secondary 
magnetospheric X-rays. Even though neither the heated polar-cap emission nor the magnetospheric emission 
are taken into account as the soft photon field illuminating the gap in this paper, they all work, in general, 
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to increase the pair creation within the gap, which cancels the initial decrease of created pairs more strongly 

than the case of the fixed soft photon field. 

Because of such negative feedback effects, solution exists in a wide parameter space. For example, 
the created current density is almost unchanged for a wide range of hm (e.g., compare the dash-dotted 
and dash- dot- dot- dotted curves in fig. 7). On these grounds, although the perturbation equations are not 
solved under appropriate boundary conditions for the perturbed quantities, we conjecture that the particle 
accelerator is electrodynamically stable, irrespective whether the X-ray field illuminating the gap is thermal 
or non-thermal origin. 

4.3. Local vs. global currents 

Let us briefly look at the relationship between the locally determined current density jgap (eq. [43]) 
and the globally required one, jgiobai- It is possible that jgiobai is constrained independently from the gap 
electrodynamics by the dissipation at large distances (Shibata 1997), which provides the electric load in the 
current circuit, or by the condition that the magnetic flux function should be continuous across the light 
cylinder, as discussed in recent force- free argument of the trans-fleld equation (Contopoulos, Kazanas, Fendt 
1999; Goodwin et al. 2004; Gruzinov 2005; Spitkovsky 2006). In either case, jgiobai will be more or less 
close to unity (i.e., typical GJ value). On the other hand, as demonstrated in § 3.3, for super-GJ cases, 
jgap is automatically regulated around 0.9 for a wide parameter range. Thus, the discrepancy between jgap 
and jgiobai is small provided that is large enough to maintain the created current density at a super-GJ 
value. The small imbalance jgiobai ^ jgap may have to be compensated by a current injection across the 
outer boundary (if the gap terminates inside of the light cylinder, charged particles could be injected from 
the outer boundary), or by an additional ionic emission from the stellar surface (if the imbalance leads to an 
additional residual at the surface). In any case, the injected current is small compared with jgap; thus, it 
will not change the electrodynamics signiflcantly, even though the gap active region may be shifted to some 
degree along the magnetic field lines, as demonstrated by Hirotani and Shibata (2001a, b; 2002), TSH04, 
and TSHC06. 



4.4. Created pairs in the inner magnetosphere 

Let us devote a little more space to examining the particle flux along the open field lines that do 
not thread the gap (i.e., hm < h < 1). Since -By vanishes, the created, secondary pairs emit synchrotron 
photons, which are capable of cascading into tertiary and higher-generation pairs by 7-7 or 7-B collisions. 
Examining the cascade, we can calculate the rate of pair creation, which takes place mainly in the inner 
magnetosphere. Denoting that the pair creation rate is K^{h) ■ {flB /2Tre) per unit area per second, we find 
Kw = 2.2 X 10*, 2.1 X 10*, 1.9 X 10*, 1.8 x 10*, 1.6 x 10*, 1.1 x 10*, 1.0 x 10*, 0.99 x 10*, 0.95 x 10*, 
0.96 X 10*, 1.0 X 10* at h = 0.20125, 0.2025, 0.20375, 0.205, 0.305, 0.405, 0.505, 0.605, 0.705, 0.805, and 
0.905, respectively, for ai ~ 70°, /i = 6.0 x 10''" G cm'', and ft-m = 0.200. Thus, the averaged creation rate 
becomes Kw = 1-4 x 10*(rii?/27re) pairs per unit area per second, giving iVpab = 3.8 x lO'^^s^^ as the pair 
creation rate in the entire magnetosphere. It should be noted that this A^pair appears less than the constraints 
that arise from consideration of magnetic dissipation in the wind zone (Kirk & Lyubarsky 2001, who derived 
10*" s^^), and of Crab Nebula's radio synchrotron emission (Arons 2004, who derived 10*°-5s-i). 

Due to strong synchrotron radiation, these inwardly created particles quickly lose energy to fall onto 
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the lowest Landau level, preserving longitudinal momentum ^ rricC per particle. These non-relativistic 
particles possess momentum flux of 2iVpair"^cc/(7^i?|(-,) ^ 1.0 x lO^^dyn cm~^, where Rpc ~ -y/rJ/njLC 
denotes the polar cap radius. On the other hand, surface X-ray field has the upward momentum flux of 
6.0 X 10^(A;T/100eV)*dyn cm~ . Thus, the created pairs will not be pushed back by resonant scatterings. 
They simply fall onto the stellar surface with non-relativistic velocities. The luminosity of e~-e+ annihilation 
line is about 3x 10^^ ergs s~^, which is negligible (e.g., compared with the 7-ray luminosity ~ lO^**-^ ergs s~^). 
On these grounds, for the Crab-pulsar, we must conclude that the present work fails to explain the injection 
rate of the wind particles, in the same way as in other outer-gap models. 

4.5. Comparison with Polar-slot gap model 

It is worth comparing the present results with the polar-slot-gap model proposed by Muslimov and 
Harding (2003; MH04a; MH04b), who obtained a quite different solution (e.g., negative E^^ in the gap) 
solving essentially the same equations under analogous boundary conditions for the same pulsar as in the 
present work. The only difference is the transfield thickness of the gap. Estimating the transfield thickness 
to be A^sG ~ hmr*\/r /wi^Ci which is a few hundred times thinner than the present work, they extended 
the solution (near the polar cap surface) that was obtained by MT92 into the higher altitudes (towards the 
light cylinder). Because of this very small A/sg, emitted 7-rays do not efficiently materialize within the gap; 
as a result, the created and returned positrons from the higher altitudes do not break down the original 
assumption of the completely-charge-separated SCLF near the stellar surface. 

To avoid the reversal of E^]^ in the gap (from negative near the star to positive in the outer magneto- 
sphere), or equivalently, to avoid the reversal of the sign of the effective charge density, = p — Pgj, along 
the field line, MH04a and MH04b assumed that Pes/B nearly vanishes and remains constant above a certain 
altitude, s = Sc, where is estimated to be within a few neutron star radii. Because of this assumption, E^^ 
is suppress at a very small value and the pair creation becomes negligible in the entire gap. In another word, 
the enhanced screening is caused not only by the proximity of two conducting boundaries, but also by the 
assumption of d{pes/ B)/ds = within the gap (see eq. [68]). To justify this p/B distribution, MH04a and 
MH04b proposed an idea that p should grow by the cross field motion of charges due to the toroidal forces, 
and that peff /B is a small constant so that cpes / B may not exceed the flux of the emitted charges from the 
star, which ensures the equipotentiality of the slot-gap boundaries (see § 2.2 of MH04a for details). 

The cross-field motion becomes important if particles gain angular momenta as they migrate out- 
wards to pick up energies which is a non-negligible fraction of the difference of the cross-field poten- 
tial between the two conducting boundaries. Denoting the fraction as e, we obtain T'mcc'^ipil{'cu/c)'^ = 
eeBAlsG (Mestel 1985; eq. [12] of MH04a). If we substitute their estimate A^sg ~ r*/20, we obtain 
e ~ 0.33(<y9/^l)77S^Ve^(w/■c^7Lc)^ where 77 = T/IO^, Be = B*/10^G, and re = r*/10km; therefore, the 
cross-field motion becomes important in the outer magnetosphere within their slot-gap model. A larger value 
of Also is incompatible with the constancy of Pes/B due to the cross-field motion in the higher altitudes. 

As for the equipotentiality of the boundaries, it seems reasonable to suppose that c\pcs\/B < c\p^\/Bf 
should be held at any altitudes in the gap, as MH04a suggested, where denotes the real charge density in 
the vicinity of the stellar surface. However, the assumption that pes/B is a small positive constant may be 
too strong, because it is only a sufficient condition of c\pef[\/B < c]p*]/S*. 

In the present paper, on the other hand, we assume that the magnetic fluxes threading the gap is 
unchanged, considering that charges freely move along the field lines on the upper (and lower) boundaries. 
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As a result, the gap becomes much thicker than MH04a,b; namely, Also ~ 0.5/iniC!7LCj which gives e < 10~ . 
Therefore, we can neglect the cross-field motion and justify the constancy of p/B in the outer region of 
the gap, where pair creation is negligible. In the inner magnetosphere, Pes/B becomes approximately a 
negative constant, owing to the discharge of the copiously created pairs. Because of this negativity of 
Pcs/B, a positive i^n is exerted. For a supcr-GJ solution, wc obtain + Jion ~ 0.9 > pcn/{^B /2Tr). which 
guarantees the equipotentiality of the boundaries. For a sub-GJ solution, a problem may occur regarding 
the equipotentiality; nevertheless, we are not interested in this kind of solutions. 

In short, whether the gap solution becomes MH04a way (with a negative as an outward extension of 
the polar-cap model) or this-work way (with a positive E^^ as an inward extension of the outer-gap model) 
entirely depends on the transfield thickness and on the resultant pcs/B variation along the field lines. If 
AZsG r*/10 holds in the outer magnetosphere, Pes/B could be a small positive constant by the cross-field 
motion of charges (without pair creation); in this case, the current is slightly sub-GJ with electron emission 
from the neutron star surface, as MH04a,b suggested. On the other hand, if Also > a7Lc/40 holds in the 
outer magnetosphere, Pcs/B takes a small negative value by the discharge of the created pairs (see fig. 5); in 
this case, the current is super-GJ with ion emission from the surface, as demonstrated in the present paper. 
Since no studies have ever successfully constrained the gap transfield thickness, there is room for further 
investigation on this issue. 
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A. Appendix 

Explicit expressions of equations (20) and (21) are as follows: 



i~i\2 /i~i\2 1 /<-\\2 

os\ ee I t'sX kq f as \ 



f ds\ ( 96** ^ . ee ( ( \ ko f ds\ ( dO^ 
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Lorentz factor 

Fig. 14. — Pure-curvature (solid) vs. synchro-curvature radiation-reaction forces. For the latter, we adopt 
the pitch angles x = 10^^, 10"^, and 10"'' rad for the dashed, dash-dotted, and dotted curves, respectively, 
and _B = 3 X 10^ and 10^ G for the thin and thick ones, respectively. Curvature radius is assumed to be 
0.4n7Lc = 6.37 x 10^ cm. 



